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Abstract 

The Schrodinger cat male and female states are discussed. The Wigner and Q-functions 
of generalized correlated light are given. Linear transformator of photon statistics is reviewed. 

1 Introduction 

The integral of motion which is quadratic in position and momentum was found for classical 
oscillator with time-dependent frequency by Ermakov fTj]. Two time-dependent integrals of motion 
which are linear forms in position and momentum for the classical and quantum oscillator with 
time-dependent frequency were found in [H; for a charge moving in varying in time uniform 
magnetic field, this was done in [Q. For the multimode nonstationary oscillatory systems, such 
new integrals of motion, both of Ermakov's type (quadratic in positions and momenta) and linear 
in position and momenta, generalizing the results of were constructed in |^. We will consider 
below the parametric oscillator using the integrals of motion. The Wigner function of multimode 
squeezed light is studied using such special functions as multivariable Hermite polynomials. 

The theory of parametric oscillator is appropriate to consider the problem of creation of photons 
from vacuum in a resonator with moving walls (with moving mirrors) which is the phenomenon 
based on the existence of Casimir forces (so-called nonstationary Casimir effect). The resonator 
with moving boundaries (moving mirrors, media with time-dependent refractive index) produces 
also effect of squeezing in the light quadratures. In the high energy physics very fast particle 
collisions may produce new types of states of boson fields (pions, for example) which are squeezed 
and correlated states studied in quantum optics but almost unknown in particle physics, both 
theoretically and experimentally. 



2 Multimode Quadratic Systems 

The generic nonstationary linear system has the Hamiltonian 

H=^QB{t)Q + Cm, (1) 

where we use 2N-vectors Q = {pi, p2,...,PN, 1i, <12i---i<1n) and C(t), as well as 2Nx2N- 
matrix -B(t), the Planck constant h = 1. This system has 2N linear integrals of motion 
which may be written in vector form 



Qo(t) = A(t)Q + A(t). 



(2) 



The real symplectic matrix A(t) is the solution to the system of equations 

A{t) = A(t)S5(t), 

A(0) = 1, (3) 

where the real antisymmetrical matrix S is 2N-dimensional analog of the Pauli matrix iCy, and 
the vector A(t) is the solution to the system of equations 

A{t) = A(t)SC(t), 

A(0) = 0. (4) 

If for time t = 0, one has the initial Wigner function of the system in the form 

W{p,q,t = 0) = Wo{Q), (5) 

the Wigner function of the system at time t is (due to the density operator is the integral of 
motion) 

l^(p,q,t) = l^o[A(t)Q + A(t)]. (6) 

This formula may be interpreted as transformation of input Wigner function into output Wigner 
function due to symplectic quadrature transform (|^). An optical linear transformator of photon 
distribution function using this output Wigner function is suggested in J^. 

The Hamiltonian (|l|) may be rewritten in terms of creation and annihilation operators 

H = ^AD{t)A + E{t)A, (7) 

where we use 2N-vectors A = ( ai, 02, . . . , a^, ctit, (12], ■ ■ ■ , ^ivt) and E(t), as well as 2Nx2N- 
matrix D{t). This system has 2N linear integrals of motion ^ which are written in vector 
form 

Ao(t) = M(t)A + N(t). (8) 
The complex matrix M{t) is the solution to the system of equations 

M{t) = M{t)aD{t), 

M(0) = 1, (9) 

where the imaginary antisymmetric matrix a is 2Nx2N-analog of the Pauli matrix — ay, and 
the vector N(t) is the solution to the system of equations 

N(t) = M{t)aE{t), 

N(0) = 0. (10) 

Analogously to the Wigner function evolution, if for time t = 0, one has the initial Q-function 
of the system in the form 

Q{a,a*,t = 0) = Qo{A), A={a,a*), (11) 



the Q-function of the system at time t is 



Q{a,a*,t) = Qo[M{t)A + N{t)]. (12) 

Here a = { q + ip) / \f2. 

For time-independent Hamiltonian (|1|), the matrix A(t) is 

A(t) =exp(E5t), (13) 

and the vector A(t) is 

A(t) = / exp(Sfir) E C(r) c/r. (14) 

JO 

For time-independent Hamiltonian (^, the matrix M(t) is 

M(t) =exp((TL>t), (15) 

and the vector N(i(;) is 

N(t) = / exp(aL)r) (T E(r) c/r. (16) 

JO 

For time-dependent hnear systems, the Wigner function of generic squeezed and correlated state 
(generalized correlated state ) has Gaussian form and it was calculated in |Q. 

Thus the evolution of the Wigner function and Q-function for systems with quadratic Hamilto- 
nians for any state is given by the following prescription. Given the Wigner function Vr(p, q, t = 0) 
for the initial moment of time t = 0. Then the Wigner function for time t is obtained by the 
replacement 

Ty(p,q,t)= W{p{t), q(t), t = 0), 

where time-dependent arguments are the linear integrals of motion of the quadratic system found 
in 0, [0, and This formula was given as integral with (5-function kernel in |]10[. The linear 



integrals of motion describe initial values of classical trajectories in the phase space of the system. 
The same ansatz is used for the Q-function. Namely, given the Q-function of the quadratic system 
(5(B, t = 0) for the initial moment of time t = 0. Then the Q-function for time t is given by 
the replacement 

Q(B, t) = Q(B(t), t = 0), 

where the 2N- vector B(t) is the integral of motion linear in the annihilation and creation op- 
erators. This ansatz follows from the statement that the density operator of the Hamiltonian 
system is the integral of motion, and its matrix elements in any basis must depend on appropriate 
integrals of motion. 



3 Multimode Mixed Correlated Light 



The most general mixed squeezed state of the N-mode light with a Gaussian density operator p 
is described by the Wigner function W^(p, q) of the generic Gaussian form. 



W{p,q) = detM"2 exp 



1 



(Q- < Q >)M^(Q- < Q >) 



(17) 



where 2N parameters < pi > and < qi >, i = 1,2, . . . , N, combined into vector < Q>, are 
average values of quadratures, 



< P > 

< q > 



Tr pp, 
Tr pq. 



A real symmetric dispersion matrix M consists of 2N^+N variances 
1 



Map = 2 < QaQp + QpQa > ' <Qa><Qf3 >, 



a,p=l,2,...,2N. 



;i8) 



(19) 



They obey uncertainty relations constraints [^]. According to previous section the Wigner function 
of parametric linear system with initial value (|1^ is 



iy(p,q, t) = detM 2 exp 

The photon distribution function of the state ([l 

Vn = Tr p|n >< n| 



i(A(t)Q + A{t)- < Q >)M-i(A(t)Q + A(t)- < Q >) 



, (20) 



n = (rai,n2, . . .,np^), 



where the state |n > is photon number state, which was calculated in [0, [|T2[ and it is 



K = Po- 



ll! 



(21) 



(22) 



The trace (^) may be calculated using the explicit form of the Wigner function of the opera- 
tor |m >< n| (see, [§) which is the product of Wigner functions of one- dimensional oscillator 
expressed in terms of Laguerre polynomials of the form 



Wran{p,q) = 2 



m— n+1 / 




n\ I q — tp 
m! \ ^2 



(23) 









-1/2 


Vo = 


det 


(m + ^l2^) 


exp 


where we introduced 


the 


matrix 





The function Hj^^(y) is multidimensional Hermite polynomial. The probability to have no 
photons is 

r / 1 XT -1/2 

< Q > (2M + l2Ar)"' < Q > , (24) 



(25) 
(26) 

(27) 



R = 2U^(1 + 2M)-^U* - aNx, 



and the matrix 





In 



The argument of Hermite polynomial is 

y = 2U*(l2Ar - 2M)-i < Q >, 



and the 2N-dimensional unitary matrix 



V2\ In 



In 



:28) 



is introduced, in which Iat is the NxN-identity matrix. Also, we use the notation 

ni\n2\ ■ ■ ■ riAf!. 



n! 



The mean photon number for j-th mode is expressed in terms of photon quadrature means 
and dispersions 



1 



(29) 



The photon distribution function for transformed state 
(IH)^(HD but with changed dispersion matrix 



It 



is given by the same formulae (p^), 

(30) 



and quadrature means 



< Q >= A-^(A- < Q >) 



(31) 



N) 



Then 



Thus we have a linear transformator of photon statistics suggested in 0. 

Let us now introduce a complex 2N-vector B = (/5i, P2, ■ ■ ■ , Pn, Pi, P2, ■ 
the Q-function is the diagonal matrix element of the density operator in coherent state basis 
I /32, ■ ■ ■ , Pn > ■ This function is the generating function for matrix elements of the density 



operator in the Fock basis |n> which has been calculated in |jT2|. In notations corresponding to 
the Wigner function (|1^) the Q-function is 



Q(B) =Poexp 



'-B{R + aN.)B + BRy 



(32) 



Thus, if the Wigner function ([T7| ) is given one has the Q-function. Also, if one has the Q-function 
(|3^), i.e., the matrix R and vector y, the Wigner function may be obtained due to relations 



M = U*(i? + ajvx)"^Ul" - 1/2, 
<Q> = V*[l-{R + aN.)-'aN.]y. 



(33) 



Multivariable Hermite polynomials describe the photon distribution function for the multimode 
mixed and pure correlated light |Tl|], |13], [1^]. The nonclassical state of light may be created due 



to nonstationary Casimir effect [|l^, ||T6[ and the multimode oscillator is the model to describe 
the behaviour of squeezed and correlated photons. 



4 Parametric Oscillator 

For the parametric oscillator with the Hamiltonian 

H -- 



+ J ■ 



2dx 



(34) 



where we take h = m = uj{0) = 1, there exists the time- dependent integral of motion found in 

i 

71 



A = ^[e{t)p - £{t)x], (35) 



where 

e{t) + iu^it)e{t) = 0, e{0) = 1, e{0) = z, (36) 
satisfying the commutation relation 

[A, A]] = l. (37) 

It is easy to show that packet solutions of the Schrodinger equation may be introduced and 
interpreted as coherent states 0, since they are eigenstates of the operator A (p5D, of the form 

^„(x, t) = ^o{x, t) exp i-^-^ + j , (38) 

where 

vl/o(x,t) = 7r-V4e(t)-V2exp!i^ (39) 

is analog of the ground state of the oscillator and a is a complex number. 

Variances of the position and momentum of the parametric oscillator in the state (PSD, (^) 

are 

(^x = ^ , (Tp = ^ , (40) 

and the correlation coefficient r of the position and momentum has the value corresponding to 
minimization of the Schrodinger uncertainty relation 

11 

If (Ta; < 1/2 (cTp < 1/2) we have squeezing in photon quadrature components. 

The analogs of orthogonal and complete system of states which are excited states of stationary 
oscillator are obtained by expansion of ( pSf ) into power series in a. We have 

*»'^-*) = (ll)'"";^*»(^-"^'»(kM)- <^2) 



and these squeezed and correlated number states are eigenstates of invariant A'^A. In case of 
periodical dependence of frequency on time the classical solution in stable regime may be taken 
in Floquet form 

e{t) = e'^'u{t), (43) 

where u{t) is a periodical function of time. Then the states (H2|) are quasienergy states realizing 
the unitary irreducible representation of time translation symmetry group of the Hamiltonian 
and the parameter k determines the quasienergy spectrum. Unstable classical solutions give 
continuous spectrum of quasienergy states. 



The partial cases of parametric oscillator are free motion 
oscillator ( u;^(t) = 1 ), and repulsive oscillator { u'^it) = —1 



[ uj{t) = ), stationary harmonic 
The solutions obtained above are 



described by the function e{t) which is equal to e{t) = 1 + it, for free particle, 
usual oscillator, and e{t) = cosht + zsinht, for repulsive oscillator. 
Another normalized solution to the Schrodinger equation 



e{t) 



for 



where 




e*(t)a^\ , y/2ax 
cosh — -— . 

e{t) ' 



(44) 



(45) 



is the even coherent state |jT8[ (the Schrodinger cat male state). The odd coherent state of the 
parametric oscillator (the Schrodinger cat female state) 



^'a/(x, t) = 2Nf'^o{x, t) exp 



\a\ 



2e{t) 



sinh 



V2, 



ax 



sit) 



where 



Nf 



exp(|a|V2) 
2A/sinh laP 



(46) 



(47) 



satisfies the Schrodinger equation and is the eigenstate of the integral of motion (as well as the 
even coherent state) with the eigenvalue a^. These states are one-mode examples of squeezed and 
correlated Schrodinger cat states constructed in |T^. The experimental creation of the Schrodinger 



cat states is discussed in ||2^. These states belong to family of nonclassical superposition states 
studied in [BTl, IE3 
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